Removing correlations in signals transmitted over a quantum memory channel 
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We consider a model of bosonic memory channel, which induces correlations among the transmitted signals. 
The application of suitable unitary transformations at encoding and decoding stages allows the complete removal 
of correlations, mapping the memory channel into a memoryless one. However, such transformations, being 
global over an arbitrary large number of bosonic modes, are not realistically implementable. We then introduce 
a family of efficiently realizable transformations which can be used to partially remove correlations among 
errors, and we quantify the reduction of the gap with memoryless channels. 
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I. INTRODUCTION 

Any physical medium that can store or transfer quantum 
degrees of freedom can be formally described as a quantum 
channel, that is, a completely positive and trace-preserving 
map on the set of quantum states. In information theory, quan- 
tum channels are mainly intended as means to convey either 
classical or quantum information [1]. They model the noisy 
interaction with the environment, causing losses, dephasing 
and decoherence. The action of the noise is usually to degrade 
the information carried by the physical system, hence reduc- 
ing the performance of quantum communication protocols. A 
standard strategy to protect the information against noise is 
the use of quantum error correcting codes |2]. These latter 
exploit redundancy and dictate the way information should 
be encoded/decoded in quantum systems. In this scenario, a 
common assumption is that the noise acts identically and inde- 
pendently at each use of the quantum channel. That is, when 
a train of signals is sent through the quantum channel, each 
signal independently experiences the noisy transformation. A 
quantum channel with this property is said to be memoryless. 

Recently, attention has been devoted to quantum memory 
channels as well (see e.g. |3] for discrete channels and |4] for 
continuous channels), where the action of the noise at differ- 
ent channel uses is either non indentical or non independent. 
From a physical point of view, that is the case when the typ- 
ical environmental relaxation times are comparable with the 
time delay between two signals. From an information theo- 
retical point of view, the presence of correlations in the noise 
can substantially reduce the efficacy of standard quantum er- 
ror correcting codes |5]. 

Then there is a twofold way to deal with information trans- 
mission over memory channels: either consider new encod- 
ing/decoding procedures accordingly to codes specifically de- 
vised for correlated errors |@], or try to remove correlations 
and then apply standard encoding/decoding procedures. Here, 
we shall investigate the second line which results unexplored 
so far We shall restrict our attention to bosonic memory chan- 
nels. 

In Ref.|0] it has been discussed a technique, called mem- 
ory unravelling, that allows the removal of correlations in a 
large class of bosonic memory channels. Actually the chan- 



nels can be traced back to memoryless ones by global unitary 
transformations prepended and appended to the channel map. 
However this is just a mathematical trick useful for evaluating 
capacites. For practical transmission rates it is unreasonable 
to realize transformations involving a large number of modes 
(channel uses). We shall hence introduce unitaries that can be 
implemented efficiently but only allow a partial reduction of 
the correlations in the memory channel. 

Specifically we shall consider the model of lossy bosonic 
memory channel introduced and characterized in 1 8] for which 
memory unravelling applies. It is known |9] that for the coiTe- 
sponding memoryless attenuating channel a standard encod- 
ing/decoding to get close to the classical capacity is given 
by (single mode) coherent states and dyne measurements 
llOI. Then, we shall evaluate the performance of such encod- 
ing/decoding procedure upon efficient, but partial, removal of 
correlations. In particular we quantify the reduction of the gap 
with memoryless channels by means of input-output mutual 
information. 

The article will proceed as follows. In Sec.|II]we shall intro- 
duce the model of bosonic memory channel under scrutiny; in 
Sec. |III]we shall consider the unitary pre-processing and post- 
processing transformations allowing the complete removal of 
correlations; in Sec.|IV]we shall introduce an efficient scheme 
for partially reducing the correlations based on unitaries act- 
ing on a finite number of bosonic modes. Numerical results 
about the reduction of the gap with memoryless channels in 
the case of chaining unitaries acting on two and three modes 
will be presented in Sec.[V] Conclusions will be drawn in Sec. 

m 



II. THE MEMORY CHANNEL 

As a case study, we consider the Gaussian memory chan- 
nel introduced and characterized in i^, which models mem- 
ory effects in linear attenuating media, acting on a set of 
bosonic modes (e.g., normal modes of the electromagnetic 
field). In the same spirit of the Refs. dH [Til], the mem- 
ory channel under scrutiny is defined as a sequence of ele- 
mentary transformations, which are concatenated through a 
memory mode. An elementary transformation involves a pair 
of ingoing bosonic modes, described through the ladder op- 
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erators {a, a'} and {c, c'}, and the corresponding outgoing 
modes {b,b''}, {d,d^}. The modes are coupled via a uni- 
tary interaction described by an operator of the form Ubs — 
exp (da) c — h.c), modeling an exchange interaction mixing 
the two modes. In the following we assume the parameter 
B to be real and positive, yielding the following Heisenberg- 
picture transformations on the ingoing modes 



^ — C^BS^ oi, [/bs — COS Q a — sin Q c , 
d = f^BS^ c C/bs — cos 6* c + sin Q a , 



(la) 
(lb) 



together with their hermitian conjugates. The quantity cos^ 9 
is usually referred to as the transmissivity parameter 

To define the action of the quantum memory channel, we 
consider a sequence of n consecutive channel uses. Such a 
sequence is associated with a collection of n bosonic modes 
with ladder operators {oj, a/}j=i....„, representing the chan- 
nel inputs, and the corresponding channel outputs associated 
with the ladder operators {foj, ...„. We model the 

channel environment as a collection of environmental modes 
with {cj, e/}j=i_...„. The memory modes, which account for 
the flux of information from one channel use to the following, 
are associated to the ladder operators {m^ , m/}j=i,...„. 
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FIG. 1: Left: a single use of the attenuating memory channel is de- 
scribed as an elementary transformation which is the concatenation 
of two beam-splitters, respectively characterized by the transmissiv- 
ities e and ry. The first beam-splitter couples the memory mode with 
the local environment, the second one mixes the memory mode with 
the input mode. Right: n uses of the memory channel are described 
as the concatenation of the elementary transformation l3L lll[l . The 
concatenation is obtained by identifying, for any fc, the outgoing 
memory mode at the fc-th channel use with ingoing memory mode 
at the (fc -I- l)-th. The memory less limit is obtained by cutting the 
information flow through the memory mode, i.e., setting e = 0. 

A single use of the memory channel is modeled as the com- 
position of two beam-splitters with transmissivities e, 77 (trans- 
missivities are within the interval [0, 1]), as depicted in Fig.[T] 
In the Heisenberg picture, it transforms the ladder operators 
according to 



erjni 



+ \Jl- ria + e)r]f 



(2a) 



& = 0? a - V(l - -V)e- -v)m, (2b) 

together with the hermitian conjugate relations. 



The action of the memory channel upon n uses is obtained 
by identifying the outgoing memory mode at each channel use 
with the ingoing memory mode at the following one (see Fig. 
[Til. In the Heisenberg picture, by iteration of Eq.s dU, we have 



fe-i 



k-j-i 



fc 



(3) 



for k — 1, . . .n. From Eq. ^ it follows that, for each k, the 
output mode operator is a function of the corresponding 
input mode operator and of the input mode operators in its 
past: Gj, with k > j. Hence, a certain amount of information 
from the j-th input flows to the fc-th output, for k > j. That 
is, the input signals from different channel uses interfere at 
the channel outputs, leading to memory effects in the quantum 
channel. 

Equation (O can be concisely rewritten in the following 
form: 



6fc = ^ fkjaj + ^ gkjCj + tkuii 



(4) 



Finally, to exhaustively define the channel model, we have 
to fix the initial state of the local environmental modes 
{cj, e|}. Different choices for the environment states lead to 
channels with different features. In the following we assume 
the local environments to be in the vacuum state. With this 
choice for the environmental states the quantum and classical 
capacities of the memory channel can be computed exactly 



The transmissivity can be related to the ratio between the 
time delay A< between to successive channel uses and the typ- 
ical relaxation time Trei of the channel environment fE]: for 
instance we may identify e ~ exp {—At/rrei)- In particular, 
the model reduces to a memoryless attenuating channel |^ 
for e = (the input Uj only influences the output bj), and 
to a channel with perfect memory ifTH] for e — 1 (all a/s 
interacts only with the memory mode). These two limiting 
settings respectively correspond to the regime At ^ Trei, and 
At <C Trei- Intermediate configurations are associated with 
values e e (0, 1) and correspond to inter-symbol interference 
channels, for which the previous input states affect the action 
of the channel on the current input ||13I1 . 



III. REMOVING CORRELATIONS 

As it has been shown in |8], a unitary transformation ex- 
ists which completely remove the correlations in the memory 
channel. Such optimal transformation can be constructed by 
considering the coefficients fkj appearing in Eq. (IH, which 
express the linear relations between the input and output field 
operators. For any n, these coefficients define an x n matrix 
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/, which admits a singular value decomposition: 



(5) 



V 



are n x n unitary matrices and A = 
ry„) is diagonal with non-negative entries. The 



where U, 

diag(77i,r/2 

unitary matrices play the role of unraveling the correlations 
In the Heisenberg picture, it maps the input and output opera- 
tors in 



(6) 
(7) 



It follows that these collective input and output field operators 
satisfy the identities 



Bk = VmAk + ^ VkigijCj + ^ Vkit. 



imi . 



(8) 



from which it is apparent that the fc-th collective output signal 
is only influnced by the fc-th input signal, i.e., the correlations 
due to the inter-symbol interference have been completely re- 
moved. Another source of correlated noise can be represented 
by the noisy terms associated to the environmental modes with 
operators ej, mi. However, if these modes are not populated 
(or, if they are in a thermal state) their presence does not in- 
duce additional correlation terms. In fact, Eq. (O can be writ- 
ten in the following form [8]: 



Bk = ^/m^k + - VkEk , 



(9) 



where {-Efc, E'^} are effective canonical field operators associ- 
ated to the environment and memory. Equation (|9|l describes 
the bosonic memory channel as the direct product of uncorre- 
cted attenuating channels, characterized by the transmissiv- 
ities {r]k}, which are in turn obtained as the singular value 
of the matrix /. In this way, we have shown that the mem- 
ory channel is unitary equivalent to a noisy channel with in- 
dependent (although non-identical) noise. This mapping has 
been used in fs*] (see also ff\) to compute the capacities of the 
memory channel. 

From an operational point of view, if the optimal transfor- 
mations are physically implemented, standard error correcting 
codes, designed for uncorrected noise, can be applied with 
high effectiveness. From this point of view, the above uni- 
taries are not only mathematical tools for computing the chan- 
nel capacities, but they are actual pre-processing and post- 
processing of information that can be used to improve the per- 
formance of standard error correcting codes. Once the uni- 
taries have to be physically implemented, one can pose the 
question of quantifying the resources needed for their imple- 
mentation. First of all, one can notice that the above unitaries 
can be constructed by combining elementary gates coupling 
at most two bosonic modes (such as beam-splitters and phase- 
shifters). According to |14], a unitary over n channel uses 
require at most n(n + 1 )/2 elementary transformations. How- 
ever, the decomposition provided by |,14J requires the action 



of beam-splitter transformations coupling arbitrary pairs of 
channel outputs. It is hence clear that to apply these trans- 
formations, one first has to wait for n channel uses and then 
start the process. For example, for the perfect memory chan- 
nel, t — \, the memory channel is unitary itself, and it can be 
perfectly inverted by waiting the n-th channel output and then 
reverse its action. 

In the following we consider a more efficient though not 
perfect procedure. Namely, we consider alternative unitary 
operators constructed by composition of elementary gates in 
a number scaling linearly with the number of channel uses. 
For any n, these pre-processing and post-processing unitaries 
only couple a fixed number I of consecutive input and output 
modes. 



IV. PARTIAL UNRAVELING OF CORRELATIONS 

We consider a scheme in which, upon n uses of the channel, 
pre-processing and post-processing unitaries are prepended 
and appended to the memory channel. We consider the case 
of canonical unitaries, hence the input modes are transformed 
as follows 



Ak ■.= ^U'^^\j a, , 

3 

Bk :=^VW,., 6,. 



(10) 
(11) 



In order to consider unitaries which can be applied in a 
efficient way, we assume that the pre-processing and post- 
processing unitaries are obtained as the concatenation of el- 
ementary unitary transformations Ui,Vt, which only couple 
t consecutive modes, in such a way that ^ — 1 of the out- 
put modes of each elementary transformation enter the next 
one. This is depicted in Figs. [2l [3] for the case of unitaries 
of ^ = 2 and f = 3. We refer to the integer I as the 
depth of the unitary transformations. Denoting /,„ the identity 
transformation acting on m consecutive modes, we then have 

(g) In-e) and 
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FIG. 2: Pre-processing and post-processing of depth 1 = 2, involv- 
ing unitary transformations which couple two consecutive bosonic 
modes at the input of the memory channel (pre-processing), and two 
consecutive output modes (post-processing). 



The application of the pre-processing and post-processing 
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FIG. 3: Pre-processing and post-processing of deptli £ — 3, involv- 
ing unitary transformations wfiich couple three consecutive bosonic 
modes at the input of the memory channel (pre-processing), and three 
consecutive output modes (post-processing). The generalization to 
unitaries of higher depth is straightforward. 



unitaries lead to the input-output relations 



Bk fkj A.J + ^ gkj ' 



tkmi 



(12) 



where fkj = Y.i,^^^^'' kifimU^^''\nj, 9kj = Hi.j^'^'^^ ki9i] 
andifc = ^;VWfc/ij. 

Then, we can optimize the choice of the unitaries of a given 
depth, in terms of a suitable performance quantifier. Here we 
would like to quantify the amount of correlations between the 
k-t\\ output mode and the /c-th input mode, and between the 
fc-th output mode and all the other input modes. In the case 
in which Z^'-^-' = U, V^^-* ~ U, we have the unraveling of 
correlations, that is, the correlations between the fc-th output 
and fc-th input are maximal, and the fc-th output is uncorrelated 
with the other inputs. Clearly, by increasing the depth of the 
unitaries, one gets higher and higher performances, since the 
case of unitaries with depth equal to I is obtained as a special 
case of the unitaries with depth equal to £' > I. 

We are going to use the mutual information as a quantifier 
of correlations. To define the mutual information, we need to 
consider a specific instance of encoding/decoding procedure 
for classical communication via the memory channel. From 
Ref. |9] we know that for the corresponding memoryless at- 
tenuating channel a standard encoding/decoding to get close 
to the classical capacity is given by (single mode) coherent 
states and dyne measurements 1 10]. 

Then, we assume the input states at the fc-th channel-use to 
be coherent states, with Gaussian distributed amplitude. Ex- 
plicitly, the input states at the fc-th use of the channel are co- 
herent states |afc), whose amplitude is chosen according to the 
Gaussian probability density distribution 



P{\ak)) ^ exp 



l«fc|' 
N 



(13) 



(Here and in what follows we neglect to explicitly write the 
proper normalization factor.) The amplitudes a^, are taken 
to be mutually independent for h ^ k. We assume the initial 
state of the memory mode and of the environmental modes to 
be the vacuum state. Under these assumptions, the fc-th output 
of the memory channel is also a coherent state, from Eq. ( fT2b . 



\Pk 



kim) 



(14) 



For extracting classical information from the quantum 
states, we consider the heterodyne measurement [1(11 de- 
scribed by the POVM elements 



\lk){lk 



(15) 



Since, in the chosen encoding scheme, the fc-th output is 1/3^), 
the probability of measuring the amplitude 7^ is 

P{^k\Pk)^Tr{E^,\Pk){Pk\) 

\{lk\l3k)? / I ^ |2\ .1^. 

= ~ exp (- 7fc - /3fc^) . (16) 

We can now compute the probability density that, given the 
coherent state \ak) is sent at the fc-th channel use, the ampli- 
tude is measured at the output. That conditional probability 
density is obtained by integrating the Eq. (fT6l) over all the ai 
with I 7^ fc, i.e.. 



P{lk\ak) 



exp(-|7fe-/3fc|2) J]p(|a,))d2c,;. (17) 

l^k 



We can also compute the joint probabiUty distribution 

P{-fk,ak)=P{lk\ak)P{W)) 

= f e^p{-\jk-pkf)llPi\ai))d^ai. (18) 

Using the explicit expression of the amplitude, (3k = 
J2i fkicti, it is straightforward to compute the Gaussian in- 
tegrals in Eqs. ( [TtI i. (fTSl l. 

From the probability distribution P{'ikiOik), it is easy to 
calculate the mutual information of variables 7^ and ak- The 
mutual information is (with a bit of abuse of notation) 

Ik - H[P{^k)] + H[P{ak)] - i?[P(7fc, ttfe)], (19) 
where H denotes the Shannon entropy, and 



P{lk) = / P{lk,ak)P{\ak))d^ak. 



(20) 



Analogously, one can compute the mutual information be- 
tween the measured amplitude at fc-th output and the ampli- 
tude of the coherent states ah with /i 7^ fc. In order to do that, 
let us denote o?^ the complex vector whose entries are the am- 
phtudes {ai}i^k- Hence, the conditional probability that the 
output measured amplitude has value a is 

P(7fc|a^)~ y"exp(-|7fe-/?fep)P(|afe))d'afe, (21) 

from which 

P{lk.av) = P{lk\dv)\{Pi\^i))- (22) 

l^k 

These Gaussian probability density distributions can explic- 
itly written in terms of the matrix coefficients fki in Eq. ( fT2b . 
The mutual information between the measured amplitude 7^ 
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at the k-th channel output and the all the input amplitudes ai, 
with I ^ k,is 

4 = H[P{jk)] + H[P{d^)] - H[P(7fe, d^)] . (23) 
Exphcitly, we have 



^'k — t; 



log2 Mil + log2 M22 - log2 det 



log2 Mil + log2 M33 - log2 det 



Mil M12 

M21 M22 

Mil Mi3 

M31 M33 



(24) 



where Af, , are the elements of the matrix 



with 



M = I Nfkk N 

Nn ON 




(25) 



(26) 



(27) 



By increasing the value of k, the mutual informations 
quickly converge to the limiting functions 



lim Ik ■ 



r := lim H 



(28) 
(29) 



Our aim is to find optimal pre-processing and post-processing 
unitaries, for any given £, that can maximize / and minimize 




FIG. 4: Left: a depth-2 canonical unitary realized as beam-splitter 
transformation. Right: a depth-3 canonical unitary realized as a net- 
work of three beam-splitters, according to Euler decomposition. 



V. RESULTS 



with 9u and 9y respectively. 

The optimization is then performed over the two angles On, 
0\!, which are the parameters of the beam-splitter trasforma- 
tions underlying the canonical unitaries, see Fig.|4l 

For the case of depth-3 unitaries, we consider canonical 
unitaries U^^^ and V'"^-* acting on three modes, parameterized 
according to Euler decomposition. 





cos 7 sin 7 
— sin 7 cos 7 

1 

COS (j) sin (, 
— sin (h cos ( 



cos 9 sin 9 
— sin 9 cos 9 
1 



(31) 



with (7,6', 4))u and (7, 6*, 0)v respectively. 

These transformations represent a small interferomenter 
composed by a network of three beam-splitters, as depicted 
in Fig. |4] The mutual informations are then optimized over 
the six angles (7, 9, (f))u and (7, 9, 0)v. 

For a given value of N, the maximum of / and the min- 
imum of /' over the beam-splitter parameters are plotted in 
Figs. |5] m as function of e and for different values of -q. It 
is worth remarking that the optimal choice of the parameters 
jointly maximizes / and minimizes /'. The figures show, as 
it has been remarked above, that better perfomances are ob- 
tained for higher value of £. 

Several comments are in order. First of all, the mutual in- 
formation enhancement and reduction are less pronounced for 
77 close to 1: as the channel gets closer to the ideal one, mem- 
ory effects tends to disappear. Secondly, if 77 is close to 0, the 
major increase and decrease in the mutual informations are 
for i — 2, while the differences in the mutual informations 
between the case with i = 2 and I — i are relatively small. 
Finally, for intermediate valure of 77, higher values of £ are 
more effective for higher values of e. 
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As examples, we present results for the optmization in the 
case of unitaries of depth 1 = 2 and £ = 3. For the case of 
depth-2 unitaries, we assume W*^^^ and V^^^ of the form 

\^ — sm t' cos 9 J 



FIG. 5: The mutual information / vs e, for N = 1. From bottom to 
top, without pre-processing and post-processing, with pre-processing 
and post-processing of depth 1 = 2, and with pre-processing and 
post-processing of depth £ = 3. From left to right, the three plots 
refer to 77 = 0.1, 0.5, 0.9. 
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FIG. 6: The mutual information /' vs e, for A'^ = 1. From top 
to bottom, without pre-processing and post-processing, with pre- 
processing and post-processing of depth 1 = 2, and with pre- 
processing and post-processing of depth £ = 3. From left to right, 
the three plots refer to 77 = 0.1, 0.5, 0.9. 

VI. CONCLUSION 

We have considered a model of quantum communication 
channel with memory, i.e., characterized by non i.i.d. noise. 
We have studied the possibility of using standard encod- 
ing/decoding strategies devised for i.i.d. errors by remov- 
ing correlations introduced by the quantum memory chan- 
nel. For the considered model, it is possible to identify pre- 
processing and post-processing unitary transformations which 
allow the complete removal of correlations, mapping the cor- 
related noises into independent (although non identical) ones. 
However, since these unitaries act globally on a train of trans- 
mitted signals, their implementation cannot be efficient, since 



one has to wait the transmission of long strings of signals be- 
fore the decoding process can start. We have hence considered 
some examples of unitary transformations with finite depth, 
involving two or three consecutive signals, permitting the effi- 
cient although partial reduction of the correlations introduced 
by the memory channel. According to suitable a quantifier 
such as the mutual information, the correlation reduction in- 
creases with increasing depth of the unitary transformations, 
but already at small depths it appears significant. 

By increasing the depth of the pre-processing and post- 
processing unitaries the reduction of correlations gets im- 
proved, as these unitaries become closer and closer to those 
allowing the complete removal of correlations. Hence, the 
depth could be used to bound the error probability of standard 
codes. 

Finally it is worth noticing that the introduced unitaries pro- 
vide an online pre-processing and post-processing procedure 
that combined with standard codes resembles the convolu- 
tional codes Iil5i1 . A deeper investigation of this parallelism 
is left for future work. 
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